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ABSTRACT 

We report results from a Chandra study of the central regions of the nearby, X-ray 
bright, Ophiuchus Cluster (z = 0.03), the second-brightest cluster in the sky. Our 
study reveals a dramatic, close-up view of the stripping and potential destruction of a 
cool core within a rich cluster. The X-ray emission from the Ophiuchus Cluster core 
exhibits a comet-like morphology extending to the north, driven by merging activ- 
ity, indicative of ram-pressure stripping caused by rapid motion through the ambient 
cluster gas. A cold front at the southern edge implies a velocity of 1000 ± 200 km 
s" 1 (M^0.6). The X-ray emission from the cluster core is sharply peaked. As previ- 
ously noted, the peak is offset by 4 arcsec (~ 2 kpc) from the optical center of the 
associated cD galaxy, indicating that ram pressure has slowed the core, allowing the 
relatively collisionless stars and dark matter to carry on ahead. The cluster exhibits 
the strongest central temperature gradient of any massive cluster observed to date: 
the temperature rises from 0.7 keV within 1 kpc of the brightness peak, to 10 keV by 
30 kpc. A strong metallicity gradient is also observed within the same region. This 
supports a picture in which the outer parts of the cool core have been stripped by 
ram-pressure due to its rapid motion. The cooling time of the innermost gas is very 
short, ~ 5 x 10 7 yrs. Within the central 10 kpc radius, multiple small-scale fronts 
and a complex thermodynamic structure are observed, indicating significant motions. 
These may be excited by the separation of the X-ray peak from the cD galaxy and 
its associated dark matter potential. Beyond the central 50 kpc, and out to a radius 
~ 150 kpc, the cluster appears relatively isothermal and has near constant metallicity. 
The exception is a large, coherent ridge of enhanced metallicity observed to trail the 
cool core, and which is likely to have been stripped from it. 

Key words: X-rays: galaxies: clusters - galaxies: clusters: individual: Ophiuchus - 
radiation mechanisms: non-thermal - intergalactic medium 



1 INTRODUCTION 

The Ophiuchus Cluster (z — 0.03) is the second-brightest 
cluster of galaxies in the 2 — 10 keV sky (Edge et al. 1990; 
see also Ebeling et al. 2002). This extremely massive, nearby 
system lies behind the plane of our Galaxy and is highly ob- 
scured, with a line-of-sight equivalent hydrogen column den- 
sity Nh ~ 3 x 10 21 atom cm -2 . Using data from the Einstein 
Observatory and EXOSAT, Arnaud et al. (1987) showed 
that the cluster has a high mean temperature, kT = 9.4l*'| 
keV (90 per cent confidence errors) and hosts a sharply- 
peaked cooling core. These findings were confirmed by later 
studies with ASCA (Matsuzawa et al. 1996, Watanabe et al. 
2001) and Chandra (Sun et al. 2007). Relatively cool, X-ray 
bright spectral components are associated with the cluster 



center, with a minimum measured temperature below 1 keV 
(Perez- Torres et al. 2009). The peak of the X-ray emission 
from the cluster is also offset from the X-ray centroid at 
larger scales (Arnaud et al. 1987). 

Ascasibar & Markevitch (2006) noted the presence of 
three surface brightness discontinuities in the central re- 
gions of the Ophiuchus Cluster (see also Fig. [T]), suggest- 
ing significant motion of the cool, dense core. Such fronts 
are common in clusters and are typically caused by merger 
events (e.g. Vikhlinin, Markevitch & Murray 2001; Reiprich 
et al. 2001). However, many relatively relaxed clusters, such 
as Abell 1795 and 2029 (Fabian et al. 2001; Markevitch, 
Vikhlinin & Mazzotta 2001; Clarke et al. 2004), also exhibit 
surface brightness discontinuities near their X-ray surface 
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Figure 1. Background subtracted, flat-ficldcd X-ray surface brightness maps in the 0.6 — 7.0 keV band, smoothed with a 2 arcsec 
Gaussian filter. The three surface brightness discontinuities noted by Ascasibar & Markevitch (2006) arc marked with arrows, (a) The 
central 8x8 arcmin 2 region covered by ACIS chip 7. Note the comet-like extension to the north/northwest and the cold front observed 
~40 kpc to the south of the X-ray peak (see also Ascasibar & Markevitch 2006). (b) The central 1.6 X 1.6 arcmin 2 region of the cluster. 
Two additional surface brightness discontinuities are observed. One is 10 arcsec (5.5 kpc) to the northeast of the X-ray peak, the other 
is 3 arcsec (1.5 kpc) to the west (see also Ascasibar & Markevitch 2006). 



brightness peaks. Simulations suggest that 'sloshing' of the 
gas even in cool core clusters can lead to visible fronts up to 
several Gyrs after a relatively minor merger event (see Titt- 
ley & Henriksen 2005; Ascasibar & Markevitch 2006; see also 
Markevitch & Vikhlinin 2007 and references within). A cool 
core in motion within the ambient ICM will also experience 
ram pressure. The impact of this ram pressure is of particu- 
lar interest for the survivability of cool cores following major 
mergers. 

Due to its location and large angular size, optical stud- 
ies of the Ophiuchus Cluster have proved challenging. How- 
ever, detailed spectroscopic studies suggest that the system 
is very massive, with a velocity dispersion a v — 1050 ± 50 
km s _1 (Wakamatsu et al. 2005), consistent with the mea- 
sured X-ray temperature and luminosity (Edge et al. 1990). 
Clear substructure in the galaxy velocity histogram is ob- 
served, indicating that the cluster is merging (Wakamatsu et 
al. 2005). The dominant cluster galaxy exhibits a substantial 
peculiar velocity, v p =-650 kms -1 , with respect to the clus- 
ter mean, cz — 9045 ± 30 kms -1 (Wakamatsu et al. 2005; 
Hasegawa et al. 2000). The presence of many nearby, smaller 
clusters and groups argues that Ophiuchus lies at the center 
of a supercluster (Wakamatsu et al. 2005; see also Hasegawa 
et al. 2000). 

Govoni et al. (2009; see also Murgia et al. 2009) re- 
port the discovery of a low surface brightness, diffuse ra- 
dio mini-halo surrounding the dominant cluster galaxy, 
2MASXJ17122774-2322108. The central cD galaxy is also 
detected at 1.4 GHz with a flux density of 29.3 mjy in the 
NVSS survey (Condon et al. 1998; see also Section ^. 1|) with 



a 45x45 arcsec 2 (FWHM) synthesized beam. There are a few 
other strong radio sources at the periphery of the cluster (see 
Johnston et al. 1981). However, their association with the 
Ophiuchus Cluster is uncertain. A variety of authors have 
searched for non-thermal X-ray emission from the cluster in 
the hard X-ray band (see Nevalainen et al. 2004, 2009; Eck- 
ert et al. 2008; see also Dunn & Fabian 2006), although the 
results to date remain inconclusive (see Fujita et al. 2008; 
Ajello et al. 2009). 

In this paper, we report detailed Chandra observations 
of the core of Ophiuchus Cluster. Section 2 discusses the 
data reduction and analysis method. Section 3 describes the 
main results and Section 4 provides physical interpretation. 
At a redshift of z — 0.03, one arcsec corresponds to 0.60 kpc, 
assuming a ACDM cosmology of Ho = 70 km s Mpc , 
fi M = 0.3, and ft A = 0.7. 



2 OBSERVATIONS AND ANALYSIS 
2.1 X-ray data reduction 

The Chandra observation of the Ophiuchus Cluster (obser- 
vation ID 3200) was carried out using the Advanced CCD 
Imaging Spectrometer (ACIS) in October 2002. The stan- 
dard level- 1 event lists produced by the Chandra pipeline 
were reprocessed using the CI AO (version 4.1.2) software 
package, including the appropriate gain maps and calibra- 
tion products (CALDB version 4.1.3). Bad pixels were re- 
moved and standard grade selections applied. The data were 
cleaned to remove periods of anomalously high background, 
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Figure 2. Co-added H, J, and K band 2MASS image of the inner 
regions of the Ophiuchus Cluster, showing the same region as Fig. 
lb. The green circle marks the position of the X-ray brightness 
peak. The radius of the circle (1 arcsec) represents a conservative 
estimate of the astrometric uncertainties. As reported by Perez- 
Torres et al. (2009), the X-ray peak is offset by ~2 kpc (4 arcsec) 
to the north of the central galaxy. 



Figure 3. 1.4 GHz NVSS image (Condon et al. 1998) of the 
Ophiuchus Cluster spanning the same region as figure[T^,. A bright 
central radio source is coincident with the location of the central 
galaxy. The spatial resolution of the radio data is approximately 
45 X 45 arcsec 2 (FWHM; see also Perez- Torres et al. 2009). The 
radio emission extends to the northwest, similar to the comet-like 
X-ray morphology. 



using the standard energy ranges and binning methods rec- 
ommended by the Chandra X-ray Center. The resulting net, 
clean exposure time is 49.5 ks. Separate photon-weighted re- 
sponse matrices and effective area files were constructed for 
each region analyzed. 



2.2 Spatially-resolved X-ray spectroscopy 

2.2.1 Spatial binning 

The individual regions for the spectral analysis were deter- 
mined using the contour binning method of Sanders (2006), 
which groups neighboring pixels of similar surface brightness 
until a desired signal-to-noise threshold is met. For data of 
the quality discussed here, the regions are small enough that 
the X-ray emission from each can be approximated usefully 
by a single temperature plasma model. We have focused our 
analysis on data from the central back-illuminated ACIS-S 
chip (chip 7). 

We initially carried out thermodynamic mapping with 
~10,000 net counts per region, to examine the thermody- 
namic structure of the hottest gas in the cluster on large 
scales. This allows us to measure the metallicity and tem- 
perature of gas with kT~ 10 keV to an accuracy of ~10 
per cent. Higher-resolution thermodynamic mapping of the 
central 40 kpc was also performed, with ~3,000 net counts 
per region, to examine the thermodynamic structure of the 
cooler cluster core. This allows us to measure the metallicity 
of this gas to better than 20 per cent, and the temperature of 
regions with kT< 5 keV to ~5 per cent accuracy or better. 
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Figure 4. Azimuthally-avcraged surface brightness profile in the 
0.6 — 7.0 keV band, centered on the X-ray peak. The profile is 
very sharp and no clear central flattening is observed. 
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In order to determine radial profiles for thermodynamic 
properties, we have also extracted spectra in concentric 
circular annuli, centered on the X-ray peak position. The 
widths of the annuli range from 1 arcsec to 40 arcsec. 

2.2.2 Background modelling 

Background spectra for the appropriate detector regions 
were extracted from the blank-sky fields available from the 
Chandra X-ray Center. The blank-sky fields were processed 
in an identical manner to the Ophiuchus Cluster observa- 
tion, and were reprojected onto the same coordinates using 
the appropriate aspect solution files. Background regions 
were chosen to match the same extraction regions as the 
science data. The background spectra were normalized by 
the ratio of the observed and blank-sky count rates in the 
9.5 — 12 keV band (due to the very small effective area of 
the Chandra telescope at hard X-ray energies, there is no 
significant cluster emission in the 9.5 — 12 keV band). The 
statistical uncertainty in the observed 9.5 — 12 keV flux is 
less than 5 per cent. 

To enable further refinement of the background mod- 
els, we have also extracted spectra from relatively source- 
free regions of the other back-illuminated ACIS CCD (chip 
5). Although cluster emission is still clearly present on this 
offset chip, it is significantly fainter than on the central 
chip (Chip 7) and the data from it can be compared to 
the blank-sky fields to identify any additional, strong fore- 
ground/background emission components. Such components 
are commonly detected in the spectra of sources lying be- 
hind the Galactic Plane. We indeed detect excess soft emis- 
sion with respect to the blank sky fields in the Ophiuchus 
data, which is modelled in the analysis. We note, however, 
that our results on the bright cluster core are relatively in- 
sensitive to the details of the background modelling. 



2.2.3 Spectral analysis 

The spectra have been analyzed using XSPEC (version 12.5; 
Arnaud 1996), the MEKAL plasma emission code (Kaastra 
& Mewe 1993), and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992). All spectral fits 
were carried out in the 0.6 — 7.0 keV energy band. The ex- 
tended C-statistic available in XSPEC was used for all fitting. 

The default spectral model applied to each spatial re- 
gion consists of an optically-thin, MEKAL thermal plasma 
model, at the redshift of the cluster. The normalization, tem- 
perature, and metallicitjQ in each spatial region are free pa- 
rameters in the fits. Since the Galactic column density de- 
termined from HI studies is high (Kalberla et al. 2005; see 
also Dickey & Lockman 1990), we include this as an addi- 
tional free parameter in the analysis, linked to vary in unison 
across all regions. We have also examined models where the 
column density is allowed to vary spatially. 

In Section 13.51 we briefly discuss results from a search 
for non-thermal X-ray emission from the cluster. In this case, 
additional power-law model components were included in 

1 Cluster abundances arc calculated with respect to Anders & 
Grevesse (1989). 



the analysis, following the method described by Million & 
Allen (2009). 

For the determination of deprojected radial profiles 
from the azimuthally-averaged spectra, we use the PROJCT 
model in XSPEC. 



3 RESULTS 

3.1 Imaging analysis 

Background subtracted, flat-fielded images in the 0.6 — 7.0 
keV band were extracted on a 0.984 x 0.984 arcsec 2 pixel 
scale. Fig. QJ, shows the X-ray image for the region spanning 
the central ACIS chip (approximately 8x8 arcmin 2 ). The 
peak of the X-ray emission lies at coordinates a(2000) = 
17 h 12 m 27.64 s , 5(2000) = -23°22 m 07.0 s . This peak is em- 
bedded in an extended X-ray bright halo, which has a comet- 
like tail extending to the north/northwest. A clear surface 
brightness discontinuity or 'cold front' is observed at a radius 
of 40 kpc to the south of the X-ray surface brightness peak 
(first noted by Ascasibar & Markevitch 2006). The overall 
X-ray morphology suggests that the cluster core is travel- 
ing approximately north to south, along a position angle of 
~ 160 degrees, with ram-pressure leading to the formation 
of the comet-like tail behind it. 

Fig. [Tp shows the central 1.6 x 1.6 arcmin 2 region of 
the cluster. Two additional sharp surface brightness fronts 
(also noted by Ascasibar & Markevitch 2006), are observed 
at radii of ~ 3 and ~ 10 arcsec to the west and northeast 
of the X-ray peak, respectively. These fronts appear to be 
sites of significant shearing and suggest significant swirling 
and sloshing within the inner, X-ray bright cluster core. 

Fig. shows the co-added H, J, K band 2MASS image 
for the same, central 1.6 x 1.6 arcmin 2 region. The dominant 
cluster galaxy is seen near the center of the field. The green 
circle in the figure marks the location of the X-ray peak; 
the radius of the circle provides a conservative estimate of 
the astrometric uncertainty in the X-ray astrometry (~ 1 
arcsec). As was also reported by Perez-Torres et al. (2009), 
the optical centroid is offset to the south of the X-ray peak 
by 4 arcsec (2 kpc). This offset presumably reflects the im- 
pact of ram pressure in slowing down the X-ray emitting gas, 
and the essentially collisionless nature of the stars (and dark 
matter) associated with the galaxy, which carry on ahead. 

Fig. shows the 1.4 GHz NVSS image (Condon et 
al. 1998; the resolution is approximately 45 x 45 arcsec 2 ; 
FWHM) for the same field of view as Fig. QJ,. A bright, 
central radio source is observed at a position approximately 
coincident with the position of the central galaxy. A low sur- 
face brightness radio feature extends towards the northwest, 
much like the comet-like trail in X-rays, although the signal- 
to-noise ratio is low. Higher-resolution and deeper multi- 
frequency radio data are required to determine the origin of 
this feature as well as the detailed relation to the complex 
central X-ray emission. 

The X-ray surface brightness profile for the cluster 
about the X-ray peak is shown in Fig. [4] The profile is very 
sharp and no clear central flattening is observed. The pro- 
file can be approximated by a power law model, SoR~ a , 
with slope a = 0.6112 ± 0.0015, or a f3 model with (3 = 
0.2685 ± 0.0002 and core radius r c < 0.5 kpc. However, nei- 
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Figure 5. Temperature map (upper left; in keV), pressure map (upper right; in keV cm -5 / 2 arcsec - 1 ), entropy map (lower left; in 
keV cm 5 / 3 arcsec 2 / 3 ), and metallicity map (lower right; in solar) for the field of view covered by ACIS Chip 7 (see Fig.[T^,). Individual 
regions follow lines of constant surface brightness, while maintaining a constant signal to noise threshold (~10,000 net counts per region). 
This results in statistical uncertainties on the temperature and metallicity of ~10 per cent or better. The electron number density is 
approximated as yj K/A, where K is the MEKAL normalization and A is the projected area. 

ther simple model provides a formally acceptable fit to the 3.2 Thermodynamic mapping 
data. 

3.2.1 Large scale properties 

Fig. \5\ shows the large scale thermodynamic structure of the 
cluster, on the same spatial scale as Fig.QJ,. From top to bot- 
tom, left to right, we show temperature, pressure, entropy 
and metallicity maps. These quantities show significant mor- 
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phological correspondence with the X-ray surface brightness 
image. The X-ray bright core contains relatively cool, low en- 
tropy gas. This gas has a high metallicity. Beyond the central 
30kpc radius, the temperature map is relatively isothermal. 

Beyond the inner 30 kpc region, no clear gradients in 
metallicity are observed. However, a large, metal-rich ridge 
spans position angles of -50 to 45 degree, at a radius of ~100 
kpc to the north of the cluster core. This ridge consists of 5 
independent, coherent, metal rich regions. The mean metal- 
licity of the 5 regions is 0.55 ± 0.04 solar, which is approxi- 
mately 1.8 times greater than the mean value of 0.30 ± 0.03 
solar determined for the 9 adjacent regions. A second ridge 
of enhanced metallicity appears to partially encircle the clus- 
ter core at a radius of ~ 40 kpc spanning position angles of 
-55 to 270 degrees at a radius <-^40 kpc. To the south, where 
the metallicity of the ridge is highest, it is approximately 
coincident with the location of the southern cold front. This 
southern ring-like ridge contains 10 coherent, independent 
regions of measured high metallicity. The mean metallicity 
of these 10 regions is 0.51±0.04 solar, which is approximately 
1.5 times greater than the mean value of 0.34 ± 0.04 solar 
determined from the 7 adjacent regions at larger radius. 

In the above analysis, the absorption was linked to vary 
in unison across all regions. We have also investigated other 
absorption modelsQ When allowing the absorption to be in- 
dependently free in each region, we obtain the map shown in 
Fig. [6] A trend towards higher absorption is observed in the 
cluster center, similar to Abell 478 (Allen et al. 1993; Sun 
et al. 2003). The innermost bin exhibits somewhat lower ap- 
parent absorption, although this may be an artifact due to 
the complex thermal structure of the X-ray emission from 
this region. A coherent structure of excess absorption ap- 
pears to cross through the cluster, with the western side of 
the field being more absorbed than the east. 



3.2.2 Properties of the inner cluster core 

Fig. [7] shows the temperature, pressure, entropy and metal- 
licity maps for the inner ~ 40 kpc radius region. The maps 
reveal steeply declining temperature and entropy as one 
moves inward toward the X-ray peak. The morphology of 
the lowest temperature gas traces a swirl structure also seen 
in the surface brightness map (Fig[lp). Interestingly, the 
X-ray brightness peak does not coincide with the pressure 
maximum: the latter is located to the north of the brightness 
peak, and trails the inferred southward motion. 

The central region of the cluster shows an enormous 
range of temperature, from less than 1 keV to greater than 
10 keV. This is only the second hot, massive cluster to show 
sub 1.0 keV gas in spectral maps (see also the study of Abell 
2204 by Sanders et al. 2009). If the cluster were at signifi- 
cantly higher redshift, this coolest gas would not be spatially 
resolved. 



The qualitative structures of the thermodynamic maps are 
unchanged using different absorption models. We note, however, 
that the exact value of the temperature and quantities derived 
thereof vary by up to 10 per cent, depending on the exact value 
of the column density. 
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Figure 6. Map of the measured equivalent, line-of-sight absorb- 
ing hydrogen column density (in units of 10 20 atom cm -2 ) for the 
region covered by ACIS chip 7 (same region as Fig. [TJi) . Other 
details as for Fig. [5] The statistical uncertainties in the measure- 
ments are at the ~ 5 per cent level. 

3.3 Azimuthally-averaged, radial profiles 

Fig. [8] shows the azimuthally-averaged, projected tempera- 
ture and metallicity profiles for the cluster, centered on the 
X-ray peak. The profiles reveal a steep central temperature 
gradient within the central 30 kpc radius. Beyond r ~ 50 
kpc, the cluster appears remarkably isothermal. The same 
central 30 kpc also exhibits a complex metallicity profile, 
which first rises steeply with increasing radius, from 0.1 so- 
lar within 0.5 kpc to almost twice solar at 3 — 5 kpc. Between 
5 < r < 50 kpc, the metallicity declines steadily before be- 
coming approximately constant for r > 50 kpc, out to the 
edge of the field. 

Fig-IUshows the deprojected temperature, density, pres- 
sure, entropy, cooling time and metallicity profiles, deter- 
mined under the assumption of spherical symmetry. We cau- 
tion that departures from spherical symmetry will have some 
effect on the electron density, pressure, entropy, and cooling 
time profiles. The deprojected profiles confirm the presence 
of a strong, central temperature gradient, which rises with 
increasing radius by over a factor of 10 within the central 
30 kpc. The minimum temperature, measured within the 
central 1 kpc, is ~0.7 keV, consistent with that of large, el- 
liptical galaxies. Indeed, it may be that this coolest gas was 
originally the interstellar medium of the dominant cluster 
galaxy, which has been stripped by ram pressure. By 30 kpc, 
the temperature has risen to 10 keV. This is the strongest 
central temperature gradient observed in any cluster to date. 

The central cooling time is very short, dropping below 
4 x 10 7 yrs. The central metallicity gradient is even more 
pronounced in the deprojected profiles, with Z ~ 0.25 solar 
within 1 kpc, a peak at Z ~ 3 solar at r ~ 5 kpc, and 
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Figure 7. Temperature (upper left, in keV), pressure (upper right, in keV cm -3 ), entropy (lower left, in keV cm 2 ), and metallicity 
(lower right, solar) maps for the central ~40 kpc radius of the Ophiuchus Cluster. Regions maintain a constant signal to noise threshold 
(~3,000 net counts per region), resulting in statistical uncertainties of ~5 per cent in temperature, for kT < 5 keV. Here, we approximate 
the volume of each emission region as fj Al, where A is the projected area and / is the line of sight distance; I = y/ r'^ nax — f'^ nin - '''max 
and r m i n are the maximum and minimum distances for all points in a region, from the X-ray surface brightness peak (see Henry et al. 
2004, Mahdavi et al. 2005). The number density is calculated from the MEKAL normalization. 



dropping again at larger radii. We caution, however, that 
the low central metallicity may in part be an artifact due to 
complex temperature structure in the innermost regions (see 
Werner et al. 2008 and references within). Beyond r ~ 50 
kpc, the temperature and metallicity remain approximately 
constant at kT ~ 10 keV and Z ~ 0.35 solar, respectively, 
to the edge of the field of view. 



Within 5 kpc of the X-ray peak, the deprojected pres- 
sure profile declines. The maximum pressure is observed 5 
kpc from the center. Our thermodynamic maps (Fig. [7p) 
show that the pressure peak lies to the north of the bright- 
ness peak. Extrapolating the pressure profile for r > 5 
kpc inward to the core, we estimate that the innermost re- 
gions are under-pressured by ~ 60 per cent. This argues 
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Figure 8. Azimuthally-averaged, projected temperature (left; in keV) and Fe abundance (right; in solar units) profiles, centered on the 
X-ray peak. The temperature profile shows a strong central gradient within 30 kpc. The central abundance is low, rises to a maximum 



value Z = 1.5^q ijg by 5 kpc, before declining again to a value Z • 
mctallicity remain approximately constant. 



- 0.35 at 30 kpc. For radii r > 50 kpc, the projected temperature and 



for substantial non-thermal pressure support in these re- 
gions, which is probably dominated by gas motions, as sug- 
gested by the multiple brightness fronts and X-ray-optical 
displacement seen in Fig. [2] Turbulent motions at the level of 
~ 500 km s _1 could in principle explain the thermal pres- 
sure deficit. However, it is also possible that non-thermal 
particles contribute to the total pressure support in the core. 
Higher resolution radio data are required to explore this pos- 
sibility. 

3.4 Analysis of the southern cold front 

A cool, dense core moving through a diffuse hot, gaseous 
halo will experience ram pressure. In detail, the pressure 
ratio at the cold front leading the core is related to the core 
velocity by (Landau & Lifshitz 1959; see also Vikhlinin et 
al. 2001) 



Po 
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where M = v/c is the Mach number of the core through the 
hot diffuse gas, and 7 = 5/3 is the adiabatic index. Here, 
index 1 refers to the ambient, hot gas and index to the 
moving, cool core. 

We have applied these calculations to the southern cold 
front at r ~ 40 kpc, which has a relatively simple geometry 
and appears to lead the overall motion of the cool core. For 
this analysis we examine the deprojected pressure profile in 
a 110 degree slice, centered on the center of curvature of the 
front. Accounting for the underlying ambient pressure pro- 
file, measured in other directions, we determine a pressure 
ratio at the southern cold front of po/pi = 1.3 ± 0.2. This 



corresponds to a Mach number of 0.58 ± 0.10 and, for an 
ambient temperature kT= 10.6 ±0.3 keV, implies a velocity 
of 1000 ± 200 km s _1 . 

The significant rotational motions that appear to be 
associated with the two other fronts discussed in Section 
3. II prevent any reliable quantitative determinations of their 
velocities. 



3.5 A search for possible non-thermal-like 
emission 

Recently, a detection of a non-thermal X-ray emission com- 
ponent from the Ophiuchus Cluster was claimed by Eckert et 
al. (2008) using data from the INTEGRAL satellite. These 
results are in apparent tension with those from Suzaku and 
Swift/BAT (see Fujita et al. 2008; Ajello et al. 2009)0 In 
general, proper temperature modelling is crucial in deter- 
mining the presence, or otherwise, of non-thermal spectral 
components in the X-ray band. As shown in Fig. [HJ a steep 
temperature profile exists in the central 30 kpc (~ 1 ar- 
cmin) of the cluster and must be accounted for. Residual 
calibration uncertainties will also affect the measured tem- 
perature and the inferred non-thermal emission properties. 
For example, a change in the globally-averaged temperature 
of ~1 keV, for a kT~ 10 keV system like the Ophiuchus 
Cluster, can change the expected thermal flux above 20 keV 
by ~ 30 per cent. Such effects currently make it difficult to 
provide any definite statements on the presence or other- 



3 A re-analysis of INTEGRAL data, also using XMM-Newton 
data to constrain the thermal emission at lower energies 
(Ncvalainen et al. 2009), is consistent with the upper limits ob- 
tained by Fujita et al. (2008) and Ajello et al. (2009). 
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Figure 9. Dcprojcctcd temperature (upper left; in keV), electron density (upper middle; in cm -3 ), pressure (upper right; in keV cm -3 ), 
entropy (lower left; keV cm 2 ), cooling time (lower middle; in years), and Fe abundance (lower right; with respect to solar) profiles. 
Profiles are centered on the X-ray peak and assume spherical symmetry. Significant asymmetries in the cluster are present, although 
the results obtained from deprojecting partial angular sectors of the cluster are broadly consistent. We caution that departures from 
spherical symmetry will have some effect on the electron density, pressure, entropy, and cooling time profiles. 



wise of possible non-thermal X-ray spectral components in 
the Ophiuchus Cluster. 

Despite these difficulties, we have searched for non- 
thermal spectral components in the Chandra data following 
the procedures outlined by Million & Allen (2009) . In detail, 
we have examined the statistical significance of the improve- 
ments to the fits that can be obtained with the inclusion of 
additional power-law components in the spectral modelling. 
The normalization of the power-law component in each spec- 
tral region was a free parameter, while the photon index was 
fixed to r — 2.0. The statistical signal from all regions was 
combined to determine the final result on the presence of 
possible power-law emission. 

We find no statistical signal of non-thermal-like compo- 
nents in the Chandra data, once the complex, underlying 
temperature structure revealed by our analysis is taken into 
account. Allowing the column density in each region to be 
fit independently, we obtain a 90 per cent confidence upper 
limit on the flux of any T — 2.0 power-law component from 
the central 8x8 arcmin 2 region covered by the ACIS-S chip 7 
(excluding only the central 5 kpc region where the tempera- 
ture drops precipitously and projection effects will be large) 



of 2.2xl0~ n ergs s" 1 cm" 2 in the 0.6- 7.0 keV band0This 
upper limit is approximately 10 per cent of the measured to- 
tal flux in the same region. Our results are consistent with 
those of Perez-Torres et al. (2009), Ajello et al. (2009), and 
Fujita et al. (2008). Our upper limits are at a similar level 
to the reported detection of non-thermal flux by Nevalainen 
et al. (2009) for the central 7 arcmin radius region. We note 
that our result depends slightly upon assumptions regard- 
ing the Galactic absorption in the direction of the cluster. 
Fixing the column density for all regions to Nh = 3 x 10 21 
atom cm -2 (the average value determined for all regions), 
our upper limit is increased to 2.7 x 10~ n ergs s -1 cm -2 . 
Due to observing energy band and field of view differences, 
direct comparisons with INTEGRAL results (Eckert et al. 
2008) are challenging. Additional hard X-ray data, for ex- 
ample from the forthcoming NuSTAR satellite (Harrison et 
al. 2005) are required to rule definitively on this issue. 



4 Our upper limits are determined from the sums of the power- 
law normalizations for the 95 regions with each region having an 
MCMC chain of length 10 4 samples after correcting for burn in. 
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4 DISCUSSION 

The comet-like X-ray morphology, large-scale thermody- 
namic maps, and clear pressure jump at the location of the 
southern cold front (Figs.[TJi,[S][7|), show that the cool, dense 
core of the cluster is moving rapidly through the large-scale 
hot, diffuse cluster medium. The velocity implied by the ob- 
served pressure jump at the southern cold front is ~ 1000 
km s _1 . The peculiar velocity of the cD galaxy of ~ 650 
km s _1 with respect to the cluster mean (Wakamatsu et al. 
2005; Hasegawa et al. 2000) also suggests significant, line-of- 
sight motion. Such speeds may indicate extreme 'sloshing' of 
the cluster core (e.g. Markevitch et al. 2001; Reiprich et al. 
2001; Tittley & Henricksen 2005; Ascasibar & Markevitch 
2006; see also Markevitch & Vikhlinin 2007 and references 
within) . 

The X-ray emission associated with the cool core is very 
sharply peaked. The coolest (kT ~ 0.7 keV), lowest-entropy 
gas lies within 1 kpc of the brightness peak. This material 
has a very short cooling time, with t coo i ~ 4 x 10 7 yr, and 
plausibly represents the stripped interstellar medium of the 
dominant cluster galaxy. 

The high relative velocity of the cool core argues that 
it was displaced from the base of the cluster potential well, 
presumably by recent merger activity. This is consistent with 
the slightly offset location of the X-ray peak with respect to 
the X-ray centroid at larger scales (Arnaud et al. 1987). 
The case for a recent major merger event is supported by 
the significant substructure in the galaxy velocity histogram. 
Given the environment of the present-day cluster, which lies 
at the center of a supercluster and is surrounded by smaller 
galaxy groups and clusters, a recent merger event would not 
be surprising. 

The degree to which cool cores are disrupted or de- 
stroyed during cluster mergers is an outstanding question in 
cluster physics, with important implications for cosmolog- 
ical studies (e.g. Burns et al. 2008; Mantz et al. 2009a,b; 
Ebeling et al. 2009). The comet-like morphology of the cool 
core, and the remarkably steep central temperature gradient 
(rising from ~ 0.7 keV to 10 keV within a 30 kpc radius) 
suggest that its outer regions have already been stripped by 
ram pressure, due to its rapid motion. Our data provide a 
dramatic, close-up view of the ongoing stripping and poten- 
tial destruction of a cool core in a rich cluster and have clear 
implications for the survivability of such cool cores. It will 
be important to examine, using constrained hydrodynamical 
simulations, the evolutionary path of such a cool core as it 
continues to be stripped and its motion slows. At this point 
it is not clear how much of the original cool core will remain 
by the time the system returns to equilibrium. 

The X-morphology of the inner 10 kpc region (Fig 
lb), the separation between the X-ray and optical/near-IR 
brightness peaks (the X-ray emission trails the cD galaxy) 
and the offset between the X-ray peak and thermodynamic 
pressure maximum (the pressure maximum trails the X-ray 
brightness peak), all suggest additional, significant motions 
within the cool core. The sharp fronts noted in Fig. lb sug- 
gest that significant rotational motion is present, and that 
the coolest, X-ray brightest gas is moving within the central 
potential. Such motion is not surprising, given the separa- 
tion of X-ray peak from the stellar and, presumably, dark 



matter potentials of the cD galaxy (Fig. 2)0 The sharp 
fronts seen in Fig lb are likely to be sites of shearing in- 
stabilities, which could lead to rapid mixing of the coolest 
gas with its surrounding environment. 

The ridges of enhanced metallicity in the cluster are 
particularly interesting. These large-scale, coherent struc- 
tures are almost twice as metal rich as the surrounding gas. 
The prominent metallicity ridge 100 kpc to the north of 
the X-ray peak, and trailing in its wake, appears likely to 
have been stripped from the cool core during its motion. 
This ridge, and the partial ring-like ridge of high metallic- 
ity material at a radius of ~40 kpc, may initially have been 
driven outward by the central AGN, possibly even during 
the same outburst. The proximity of this ring to the south- 
ern cold front argues that this metal rich gas is currently 
being stripped from the cool core. A qualitatively similar 
high-metallicity ridge is also seen in the central regions of 
the Perseus Cluster (Sanders et al. 2005). 

Govoni et al. (2009) show that the cD galaxy is sur- 
rounded by a bright, diffuse mini-halo of synchrotron radio 
emission. Despite the fact that a significant non-thermal par- 
ticle population must be present to explain the radio emis- 
sion, our study provides no evidence for non-thermal X-rays 
from the central regions of the cluster. High-resolution radio 
data for the cluster core are not currently available. It will 
be important to obtain such data and examine the impact 
of the cluster dynamics on the properties of the central ra- 
dio source, and the interaction between this source and its 
surrounding environment. 
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